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Abstract 

We propose a new discrete approximation to the Einstein equations, based on 
the Capovilla-Dell-Jacobson form of the action for the Ashtekar variables. This 
formulation is analogous to the Regge calculus in that the curvature has sup- 
port on sets of measure zero. Both a Lagrangian and Hamiltonian formulation 
are proposed and we report partial results about the constraint algebra of the 
Hamiltonian formulation. We find that the discrete versions of the diffeomor- 
phism constraints do not commute with each other or with the Hamiltonian 
constraint. 
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1 Introduction 



Einstein's equations are beautiful because they capture the phenomena of grav- 
itation in a simple geometrical statement: the vanishing of the Ricci curvature 
in empty regions of spacetime. However, when we use them to solve a practical 
problem in gravitation we discover they have another side — they are compli- 
cated nonlinear partial differential equations. In special cases, most typically 
when symmetries have been imposed, one sometimes can make use of the ge- 
ometry to help discover the solution to a physical problem. But, when one is 
studying the generic problem of constructing solutions to, or evolving, the Ein- 
stein equations, little of the geometric beauty comes out in the techniques we 
use to try to solve the theory. 

This is especially the case for numerical approximation methods. Such meth- 
ods are crucial for making progress with important astrophysical questions such 
as gravitational wave production by realistic sources. By the time all of the 
elements necessary for making the numerical problem well defined are in place, 
including gauge fixing and finite differencing schemes, very little of the geomet- 
rical beauty of the equations remains. 

For over thirty years there has been available an alternative to the finite 
differencing approaches to the Einstein equations, which is the Regge calculus jl]]. 
In this approach a large and generic set of solutions is constructed by limiting 
solutions to manifolds in which the curvature is restricted to have support on 
sets of measure zero. Typically, the manifold is broken up into simplices, and 
the curvature is restricted to lie on the boundaries at which simplices are joined. 

The idea of Regge calculus is that such simplicial manifolds could approxi- 
mate a given smooth manifold arbitrarily well. Unfortunately, at least up until 
this time, Regge calculus has not been developed into a powerful tool for use 
in realistic calculations. Or, at least, one should say that this is the case in the 
classical theory, as recently a version of Regge calculus has been shown to yield 
a very effective method for calculating the path integrals in quantum gravity in 
two |, |, |, |, §, three § and even four dimensions g § [K))(| 

In this paper we introduce a new discretization of Einstein equations that 
results from applying the exact Einstein's equations to a special restricted class 
of geometries. This new formulation has two basic features. First, it is based on 
the Ashtekar formulation of general relativity |f6|, [Tt], |f8| , in which the dynam- 
ical variables are frame fields and self-dual connections. Second, the restricted 
class of geometries we study are those in which all the fields are distributional. 

These are connected in the following way. As the Ashtekar formalism is 
polynomial at both the Lagrangian and Hamiltonian level, the field equations 
allow a wider class of solutions than the conventional form of the Einstein equa- 
tions. These include solutions in which the determinant of the metric vanishes, 

4 For other work on four dimensional quantum gravity using Regge calculus, see [l^ , 
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for which the usual relations for the Christoffel symbols in terms of the metric 
components would not be defined. This is possible because one of the basic 
fields of the formalism is a three dimensional frame field, which is chosen to 
have density weight one. (Actually, such solutions are allowed for all first or- 
der formulations of the Einstein's equation, as has been emphasized recently by 
Horowitz @.) 

Among such degenerate solutions are those for which the densitized frame 
field is actually distributional. Because of the Yang-Mills like gauge invariance 
of the theory, it turns out that there are many solutions in which the frame 
fields have support on one dimensional curves in the Hamiltonian formalisms. 
This is because the Gauss's law constraint is solved by frame fields which are 
covariantly divergence free. These configurations may be taken to be of the 
formQ 

E a a \x)=a 2 J ds6 3 (x,a{s))d a (s)el l (1) 

Here a is a closed or open curve in the spatial three manifold, which we will 
denote E and e 1 is an element of the Lie algebra of SU (2) which is associated 
with that curve, a is a constant with dimensions of length, which is necessary 
if, as is natural, the Lie algebra element e % a is dimensionless, so that the frame 
field can also be dimensionless. 

In the Ashtekar formalism it is always important to keep track of the density 
weights. Thus, note that the right hand side of (1) is naturally a vector density, 
as is required. 

We can consider distributional geometries of the form of (1) for complicated 
graphs or lattices. For example, let T be some graph in E, with edges 77 where 
1=1, ...N. Then we consider distributional configurations of the form 

E«(x) = a 2 £ I dsS 3 (x, 7 /(sM(s)4 (2) 
/ J 

Associated with the graph T is then a subspace of the configuration space of the 
theory which is then given by the N lie algebra elements e}. 

Now, such a distributional frame field may not seem very physical. Indeed, 
that was the first impression when expressions of this form first arose in the 
quantum theory in pi] . However, we have recently understood^], |2^, Q 
that such distributional frame fields can actually approximate arbitrarily well 
any smooth metric, q a b, as long as the scale of the curvature of the metric is large 
compared to the spacing between the links of the graph, where both distances 
are measured in terms of q a b- 

5 We use throughout index conventions of reference |E3|. Early latin indices (a, b, c, ...) are 
three dimensional spatial indices, middle latin indices (i,j,k,...) are three dimensional frame 
indices, tilde's indicate density weight one and capital middle latin indices I, J < ... refer to 
individual loops in a set of loops. Greek letters are generally used to denote loops. 
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Let us sketch here briefly the main ideas behind this, for more details the 
reader is referred to |^2|, ^3|, |24| . Certainly, many observables on the configura- 
tion space of general relativity cannot be defined on configurations of the form 
of (2). Included in these is the metric at a point. Since the metric is defined 
through the expression 

f\x) = E ai (x)E>>(x) (3) 

it is not, at least naively, well defined for configurations of the form of (2) as 
the product of the distributions is not defined. Moreover, it can be shown that 
the situation cannot be improved by any regularization or renormalization pro- 
cedure of the type that is usually used in quantum field theory. The problem is 
that any such procedure either introduces spurious dependence on the regular- 
ization procedure or changes the density character of the observable p^, [23] |. 

However, neither in classical nor in quantum physics do we ever observe 
the metric at a point of space or spacetime. Leaving aside the difficult issue 
of diffcomorphism invariant observables in general relativity, it is obvious that 
what is really observed in any real experiment in general relativity are averages 
of the metric, smeared over some regions of space or spacetime. 

The question is then whether there are observables which measure the metric 
smeared over regions of space or spacetime which can be defined on distribu- 
tional configurations of the form of (2). It turns out that there are a number of 
such observables. Among them are the following three [^2[ 

1) The area of any given two dimensional surface, iS, in S. It is denoted 
A[S). 

2) The volume of any three dimensional region 1Z in S. It is denoted V[TZ]. 

3) An observable that measures the integrated norm of any one form to on 
E. Written in terms of the classical three metric q a b it is 

Q[u]= / yj ' det{q)q ab uj a u b . (4) 

•J S 

Note that the square root is a density and is thus integrable. 

The key point is that in spite of the distributional character of the frame 
fields (2) the areas and volumes associated with them are finite when they 
are nonvanishing. For example, the area that a surface S acquires from the 
distributional configuration E^ 1 is nonvanishing as long as the surface intersects 
the graph at least once. In that case the area is 

A[S]^a 2 J2^ + [S,li]\ei\ (5) 
i 

where |e/| is the norm of the SU(2) Lie algebra and X + [5,7/] is the unoriented 
intersection number, that simply counts positively the number of intersections 
of the surface and the curve. The area observable will be described in detail in 
the appendix. 
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Thus, if we are measuring, not the metric at a point, but the areas of surfaces, 
it is possible to arrange the graph and choose the e\ such that the areas the 
surfaces have in the distributional configuration (2) approximate arbitrarily well 
the areas they have given a smooth metric q a b- For example, one way to do this 
is the following [p2|. Fix |ej| = 1. Then space the lines in the graph so that, on 
the average one line crosses every surface once per a 2 units of area of the surface, 
as it is measured by q^. This can be done consistently for any set of surfaces, 
as long as their radii of curvature are small compared to a, where again, the 
radii of curvatures are measured with respect to q a b- Then, consider any such 
surface whose area from q a i> is large compared to a 2 . Its area according to the 
distributional geometry (2), from (5), will then be equal to the value from the 
smooth metric, up to errors of order a 2 divided by that area. 

Similar statements can be made concerning the second and third observables. 
We may conclude that as long as we only measure such smeared observables, 
any smooth three metric may be approximated arbitrarily well by distributional 
metrics of the form of (2). 

Having established this correspondence between distributional and smooth 
three metrics, we may go on to ask four additional questions. First, can we de- 
fine a corresponding set of distributional connections such that the constraints of 
general relativity can be extended to the case of such distributional connections 
and frame fields? Second, can we find solutions to the constraints among these 
distributional initial data, yielding an extension of the physical phase space of 
general relativity to such distributional solutions. Third, can we define Poisson 
brackets on the space of distributional frame and connection fields that, together 
with the extention of the constraints to these fields, defines a constrained Hamil- 
tonian system? Fourth, can we find corresponding distributional solutions to 
the full set of Einstein's equations, either by evolving the distributional solutions 
to the constraints in time or by directly solving the Einstein equations. 

The purpose of this paper is to show that the answer to the first three ques- 
tions is yes. We find that the equations of general relativity can be consistently 
reduced to a set of equations that govern the evolution of distributional frame 
and connection fields. The starting point for this formulation of dynamics is 
to begin with a form of the action first written down by Capovilla, Dell and 
Jacobson|2^]. This allows us to make a consistent truncation of the constraint 
equations onto a finite, but arbitrarily large, dimensional phase space, that can 
approximate arbitrarily well smooth solutions to the constraints of general rel- 
ativity. As in full general relativity, the Hamiltonian is a linear combination of 
constraints, and those constraints are in correspondence with the constraints of 
the continuum theory so that there are discrete analogues of the Hamiltonian, 
diffeomorphism and SU (2) gauge constraints. 

The key to these results is to invent connection fields that are associated with 
the faces of the graph T. If we label the faces by an index a, the connection 
will also be specified by associating a Lie algebra a a to each face. Associated to 
every graph T we will then have a phase space Vt, in which SU(2) Lie algebra 
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elements corresponding to frame fields of the form of (2) are associated to each 
link and the conjugate fields are Lie algebra elements associated to each face. 

Given a formulation of a discrete approximation to general relativity as a 
constrained Hamiltonian system, we would like to use its evolution equations 
to construct solutions that approximate smooth solutions to general relativity. 
In order to do this, the key question which must be answered is that of the 
algebra of the constraints of the discrete theory. One's first expectation is that 
the constraint algebra is most likely second class, as a result of the fact that 
the diffcomorphism invariance of the continuum theory appears to be broken by 
the reduction to the discrete theory. This is known to occur in several attempts 
to construct Hamiltonian formulations of the Regge calculus [^6|, |28|, |29|. 
We have studied the constraint algebra in the present model and we report 
here partial results concerning its form. While we have not yet been able to 
calculate the full algebra, we have been able to compute the algebra in a certain 
limit, in which the SU(2) internal gauge symmetry is reduced to an abelian 
algebra which is £/(l) 3 . In the continuum theory this limit has been studied, 
and corresponds to the limit in which Newton's constant G is taken to zero 
in the Ashtekar formalism (3^]. It corresponds to a chirally asymmetric theory 
which includes the full self-dual sector of the theory but only the linearization 
of the antiself-dual sector. What we find is that in this limit the algebra of 
the constraints analogous to the gauge and spatial diffeomorphism constraints 
is not first class in this discrete approximation. 

Finally, it may also be interesting to construct discrete solutions to general 
relativity which approximate smooth solutions by restricting the four dimen- 
sional equations of motions of the theory to distributional fields so that time, as 
well as space, becomes discrete. We show in section 4 below that there is such a 
formulation, which is based as well on the Capovilla-Dell-Jacobson form of the 
action. 

This paper is organized as follows. In the next section we derive the dis- 
crete Hamiltonian formulation by reducing the Capovilla-Dell-Jacobson form of 
the action for general relativity to a form appropriate for fields that are distri- 
butional on three surfaces, but continuous in time. In section 3 we study the 
algebra of constraints of the resulting Hamiltonian formulation, while section 
4 is devoted to the construction of a four dimensional lagrangian in which the 
fields are distributional on the spacetime manifold. Our conclusions and sugges- 
tions for further work are in section 5. The appendix contains technical details 
about the extension of the area observable to distributional fields. 

2 The Hamiltonian formulation 

The best way to insure that a constrained Hamiltonian formulation is consistent 
is to derive it from an action principle. In this paper we will thus describe 
two closely related formulations of an action principle for distributional fields 
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in general relativity^ In the first, spacetime is discretized, so that the full 
four dimensional Einstein equations are replaced by a finite set of difference 
equations. This is the subject of section 4 below. In the second, space is 
discretized, but time is kept continuous. This is necessary if we are to derive 
a discrete approximation to the Hamiltonian theory, and is the subject of the 
present section. 

Both types of Lagrangians are constructed by using a form of the Lagrangian 
for general relativity first written down by Capovilla, Dell and Jacobson. It will 
be helpful to first explain this form, as it is the starting point for all the devel- 
opments of this paper. On spacetime, which will be denoted M, we consider 
two independent fields. These are an SU (2) connection one form, denoted A 1 
and a three by three matrix of scalar fields, denoted (jA. We may then consider 
the action 

S{A,cf) = [ F'AP'l^ (6) 



In this form of the variational principle, <jf ^ is to be varied respecting two 
constraints. These are a symmetry condition, 

<f>ij = <Pji (7) 

and a trace condition 

J2 <t>u = -6A (8) 

i 

where A is the cosmological constant (which may be zero). 

Let us assume that M. has the topology of £ x R, for some three manifold 
E, and consider a 3 + 1 decomposition of M. We may then split the spacetime 
coordinates, x a into a time coordinate t and spatial coordinates x a on S. It is 
then easy to show that if one defines the corresponding Hamiltonian theory for 
evolution in t the Ashtekar form of the Hamiltonian theory is found. 

Let us sketch this, as we will shortly be following the same procedure in 
the distributional case. We may define the frame fields E m as the conjugate 
momenta to A l a , according to the usual 

E ai = = 2U- 1 ]\B aj (9) 
SAi 

where B al = \e ahc Fl c . It is then one line to show that the action may be written 

as 



S= dt 



E ai A\ - A\Q l (10) 



6 A preliminary version of the formulation of this section was presented in [B4j|. We may 
note that several of the equations of that presentation have been corrected here, including the 
relationship between the connection and curvatures. 
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where Q l is the Gauss's law constraint defined by 

G l (x) = V a E al (x) =0 (11) 

where V a is the SU{2) gauge covariant derivative. The Poisson brackets may 
be read off from (10), they are, 

{A bj (x),E ai (y)} = 5 3 (x,y)5Z5 i j (12) 

However, because of the conditions (7) and (8) we put on faj, not all the E al (x) 
defined by (9) arc independent. Instead there are primary constraints. It is not 
hard to show that there are exactly the momentum and Hamiltonian constraints 
of the Ashtekar formalism: 

C a = F* b E H = (13) 
C = e ljk Fl b E a i E hk + Ke ljk e abc E al E b ^E ck = (14) 

where F* b is the Yang-Mills curvature associated with A % a . 

It is intriguing that in this formulation the dynamical constraints of general 
relativity are recovered as primary constraints. This will turn out to be a great 
help when constructing the discretization. 

Now, we are going to construct a discretization by restricting this form of 
the action to distributional fields. The key question to be answered is what 
kind of distributional connections and curvatures will be employed to do this. 
In order to answer this question, let us first note that by (9) the support of E al 
is likely to coincide with the support of B al = ^e abc F bc . Thus, returning for the 
moment to the three dimensional formalism, we will expect to have curvatures 
of the form, 

B a J(x) = |E/ dsS^xmismim. (15) 

Here, the G is put in for dimensions. We would like the free factors fe. to be 
dimcnsionless. In Ashtekar's formalism it is GF* b that has the dimensions of 
curvature, which is inverse length squared. 

Thus, in the Hamiltonian theory the spacetime self-dual curvature will be 
represented by one SU(2) Lie algebra element, b\, associated to each line of the 
graph. 

The question we must now ask is what form of a distributional connection 
gives such a curvature? It turns out that the right answer is that the connection 
has support on the faces of the graph. In order to make the combinatorics simple 
and explicit, for the rest of this paper we will restrict attention to cubic lattices. 
Thus, given an arbitrary coordinate chart on E let us define a standard cubic 
lattice with coordinate lattice spacing a. The vertices will be labeled by three 
integers n and the links by the pair (ft, a). Thus, "/na(s), with s G (0, 1) will be 
taken to refer to the link leaving the vertex n in the positive a direction. From 
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Figure 1: The lattice variables e\~ and aLx., where x means the direction of 
e fix respectively the normal of the face Sftzy- 



now on, a, 6, c will be one of three positive directions in a right-handed system, 
see fig(0). The vertices themselves are located at points x 7 -% — 7^ (0). 

We will then take the frame field and curvature to be of the forms (2) and 
(15), with the index / over the links now labeled by (n,a), 

Let us also label the faces by (n, ab) with a > b. Then, we may consider 
distributional connections of the form, 



3 E J ^Si^e^ix^a))^ (16) 



where c is the positive normal of the face spanned by a and b, see fig(Q), and 
a stands for an arbitrary pair of coordinates on the two dimensional surface. 
We must first compute the associated curvatures and show that they are of the 
form of (15), so that the b % n h can be expressed in terms of the a l m . To compute 
the first, derivative term of (15) smear the distribution with an appropriate one 
form fd and compute, 



E 6 ^/ ^a(s)/ai(7na(s)) (4c 
be 



(n—b)t 



(17) 

The second, non-abelian term is also found to have support only on the links 
of the lattice. The result is that the curvature of the connection (16) is of the 
form of (15), with 



E eSSe f(4a-^_ 5) a) + r ij7! (4s 



b,c^a 



+ a' 



(ah 



at 



b (n-c)bj \ nc (n-b)c 



(18) 
see also fig(|2j). 

In the continuum, the uncontracted Bianchi identity gives a relation between 
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the curvature B a% and the connection A l a , namely 

V [a (e bc]d B di ) = (19) 

It is easy to verify this equation directly by substituting 

B ai = e abc F b \ = e abc (d b Al - d c A\ + e ijk A{A k c ) (20) 

into the left-hand side of (19). Substituting the distributional fields (15) and 
(16) into the left-hand side of (19) yields 

V [a (e bc]d B di ) = £ 5 6^-6? 3 _a)a (21) 





-(b+c))a " T U (n-b)a i I V"" S + 



One might expect that substituting (18) into (21) would make the expression 
identically zero. This turns out not to be the case. The right hand side of (21), 
expressed purely in terms of a Ra via (18), has terms that are linear, quadratic, 
and cubic in at.. Only the linear and quadratic terms cancel, leaving an ex- 
pression cubic in a%~ that, in general, is not zero. 

To understand this, we remember that the Bianchi identity can be inter- 
preted geometrically as a manifestation of the "boundary of a boundary = 0" 
principle |3l|| . Taking a cube in the dual lattice with lattice vertex n at the 
center, we parallel transport a test spinor around each of its faces in such a way 
that each link gets traversed twice, once in each direction. The resulting change 
in the test spinor is, to lowest non-zero order in 6^ a , exactly given by the linear 
and quadratic (in a l R& ) part of the right hand side of (21). Thus, the Bianchi 
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identity for our distributional frame fields is 



E^-^- a)a (22) 

a 

+h ijk (4a + E [^i-(S + a))a + ^ ) + /? ( Va)a)] = 



where 



4a = E - (23) 



That is, is equal to the linear terms of b l n& in (18). 

This completes our specification of the kinematical structure of the dis- 
cretization. The next question to be addressed is how to extend the constraint 
equations to distributional fields of the form of (2) and (16). This is not a com- 
pletely trivial problem, as the computations involve products of distributions. 
For this reason, the simplest approach is to use the Lagrangian approach to 
define the dynamics. 

In the continuum theory, the action can be written 

S(A a , A a ,c/ ) ) = 2jdtJ^ [A\ - V a Al\ B^r. 1 (24) 

We now make the coefficients of the distributional fields, et. and a^-, time 
dependent and plug these fields into the above action. When we do this, we 
notice that in order to perform the integrations, some kind of restriction has 
to be put on (f>ij. To see what restriction we should choose, notice that in the 
Ashtekar formalism there are 18 configuration variables at each point in £ ; E m 
and A\ (we do not count A l Q as a configuration variable, since it plays the role 
of a Lagrange multiplier). We would therefore expect to have 18 configuration 
variables associated with each vertex of the lattice. This turns out to be the 
case if we restrict <pij (x) to be constant on the three links associated with each 
vertex. That is, 4>ij(x) has the same value, which we will denote as <pij{xf{) 1 on 
the three links obtained by starting at the vertex n and moving in the three 
positive directions. Note that (f>ij(x) is discontinous at each vertex. We now 
have 18 configuration variables associated with each vertex (9 a's and 9 </>'s). 

If we make this restriction, keeping Aq(x) unrestricted, then the discrete 
version of the action is 

,Ao,4>] (25) 



J dt£[a, 



where 



C[a, A , 4>] = 4c x {^/(zn) + %- e) c^/(^-a) (26) 
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+ 12 [^iri+a+S)A 1 ( x a+a+h) + ^+a)c^/(^+a)+ 

a,b^c a^b 

+ \ V \n-c+a+b)c^ 1 ( X n-c+a+b) + ^{n-c+a)^ (Xft-c+a)] } + 

+ ^12 A h( X n)Gk 

n 

Here, Q % n is the discrete version of the Gauss law constraint; it is given by, 

Qk =2Ea (4^ 1 (^)-^- a)a ^ 1 (^-a)) (27) 
+r ikl + &k_ a)a 0«Wa)) 



I 4a + E t^fn- 
\ b,c^a b^c 



(S+c))o + a (n-b)a] 



Following the continuum theory, let us define the coefficients of the frame 
field to be 

4a = 2 *« 1 (**)4i (28) 

The Lagrangian is now 

£[a,A>,e] = ^4><{4 + e M^ (2 9 ) 



n,c 



a,b^c a^b 
+ 2 6 (™- £ + a +&)c + e (n-c+a)c\ } + 

n 

Now, the Gauss constraint is 

Gh = 12 ~ in-a)a) + (30) 



\ b,c^a b^c 

Notice that, at each vertex, we have replaced the five independent variables 
hj{xn) with the nine variables e l na . Due to conditions (7) and (8) on faj, we 



L 2" , (™-(t+c))a + a (n-b)a) 

bjC^a b^c 
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have 4 primary constraints on et & at each vertex: 

Cfta = e a,hc e ffifiiie (31) 

Cn = ^jk^abc^na^bfic (32) 

These are the lattice vector and scalar constraints, respectively. 

In the continuum, the abelian term of the vector constraint is equal to the 
abelian term of the Gauss constraint times the connection. The physical inter- 
pretation of this is that the diffeomorphism constraint vanishes in the abelian 
limit when the G Newton constant goes to infinity |^2| . This is not the case for 
the discretized constraints (|3l"l). 



3 The constraint algebra 

We would now like to find the Poisson bracket algebra of the constraints. Notice 
that the definition of the momenta conjugate to both et~ and at. are actually 
primary constraints, since they do not involve time derivatives of el- and at a : 



dC 1 



dat. 8 



_ f t , t 

q t e nc ~i e (n—c)c 



(33) 



T \-e* - ■ 

/ j L 9 (n+a+Mc 



l 



L 2 C («+a+fc)c e (n+a)c + 



2 {n-c+a+%)c ^ (f 



-c+a)c] 1 



na 

Let us label these constraints by 

XI = (35) 

where I is an index labeling the above constraints. Let us now treat a L x M x 
iV cubic lattice, with opposite points on the boundary identified, so that the 
topology is that of a 3-torus. This means that there are (N — 1)(M — 1)(L — 1) 
independent vertices. We also require L, M and N to be even numbers. It is 
easy to see that for any particular constraint xii there exists another constraint 
Xj such that Xj} does not weakly vanish. Therefore, all constraints xi are 
second class constraints. We now follow the procedure first set forth by Dirac 
in dealing with second class constraints. We define the Dirac bracket of two 
functions on the phase space, £ and ip, as 

U, ^}db = U, ^} P b - U, Xi}pb^ IJ {x.J, ^} P b (36) 
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where fi /J is defined by 

{xi,Xj} P bCl JK = 5? (37) 

Therefore, the Dirac bracket of any function on the phase space with any \i 
vanishes by construction. 

It is not difficult to see that the resulting Dirac brackets are non-local, in 
the sense that an a 1 on one face of the lattice can have a nonvanishing Dirac 
bracket with an e l on a link arbitrarily far from it. This happens because to 
find the Dirac brackets we have to invert the matrix {xitXj}- This may seem 
unphysical, but it is actually necessary if the Dirac algebra of the a l, s and their 
conjugate momenta (defined by 33) is to be local. The problem is that the 
relation between the e*'s and the conjugate momenta are nonlocal. (Actually, 
a natural rearrangement of the configuration variables makes the Dirac bracket 
local; see below.) 

The straightforward way to proceed is to write a computer program to invert 
the relation (37) and find the Dirac bracket for any given lattice, which we did. 
However, a small trick gives us an easy form of fi. We rearrange the constraints 

= n 4."i: < 38 > 

in such a way that all new constraints contain only one et & , and the rest in the 
same manner. 

Using the above rearrangement of constraints, (37) is simple enough to solve 
analytically. The Dirac bracket is 

{<4p4J d6 = (39) 
{-l) k+l+m -\sign{k) - 5 kfi )(sign{l) + S lfi ) x 
x (sign(m) + S mfi )S a jSij 

rh — n = ka + lb + mc: a • b = a • c = b • c = 



where 



sign{l) = for 1 = 0, +1 for I > 0, -1 for l<0 (40) 



This expression was checked by directly performing the inversion required to 
solve (37) on a computer using a 6x6x6 cubic lattice. 

This expression seems nonlocal at first glance. However, if we rearrange the 
configuration variables as 

e na ~ e n-a,a e na (41) 

a na = a Ha+ ^2 i'b+&))& + a \n~~b)a ^ 

b.c^a b^c 
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Figure 3: The Dirac brackets are local for a certain combination of our original 
discrete configuration variables. 



(this will not change the number of configuration variables), we get in fact a 
local Dirac bracket, see also fig(pl): 



{a' nh ,e'rha\ db = t & :b S i,A 



(43) 



We now turn to the problem of finding the algebra of the constraints. To 
simplify the job of finding the Dirac bracket algebra of the constraints, we 
considered first an approximation in which we neglect all non-abelian terms in 
the constraints, so that 

Cna = e^bU (44) 
C H = eyfce a fcc e «a4^c (45) 
^-E( e ^-e| s _ a) a) (46) 

d 

are the vector, scalar, and Gauss constraints with 

bk = E ^ ^ - °(M>0 ( 47 ) 

As mentioned in the introduction, this approximation corresponds, in the 
continuum theory, to the limit in which Newton's constant is taken to zero, 
as described in |3(J. Note that both Bianchi identities ((TTJ) and ( |23| ) are now 
fulfilled, as there are no cubic terms at all. 

The results of the Dirac bracket algebra is now as follows: The algebra of 
a generic lattice (as specified above) does not close. Almost all dirac brackets 
between the vector and scalar constraints can't be expressed as a linear combi- 
nation of constraints. This is obvious when computing a particular bracket and 
has been checked in a linearization program in the case of a 6 x 6 x 6 lattice. 
Treated as (secondary) constraints, the non vanishing dirac brackets actually 
are first class. The algebra must eventually close, as there are a finite number 
of combinations of a^ a 's and e\ h 's. The full algebra including the brackets of the 
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new first class constraints is quite complicated for a generic lattice. Therefore, 
we present the full algebra for a simple but not trivial lattice and only a part 
for the generic lattice. 

Let us start with the simple case. For a2x2x]V lattice, meaning N — 1 
independent vertices on a line, much is simplified. Let the line be along the 
x-axes. Then, there are only N — 1 constraints left (as b Ax = which means 
that e Ax = 0), namely 

C Hx = D H = ei y bi z - 4 Z *4 for n = 1, 2...N (48) 

Now, 

{D n , Dft +mi } = (-l) m - 1 sign(m)(e z ny b t 77lz - e^tf^ + e l Ay b\ z - e^ z 6^|9) 

= ^ (n ,n-\-mx) 

, where sign(l) is defined in (f4(i|). The full algebra then becomes 

{-0(i!+li,n+mi)i ^(n+pi,ff+g£ ) } = ( — 1)' P 1 Sign(l — p)-D(H+;£ i n+px) + (50) 

+ (_iy-9-l s ^ n (j _ + (51) 

+(-l) m -P- 1 s£gn(m-p)D (f?+TO £ )S+p;3i ) + (52) 

+ (-l) m ~ q ~ 1 sign(m - q)D {n+mS:fi+qS . ) (53) 

where 

2D A = D (Htfi) (54) 
Among the (y) number of -D(n. m )'s there are only one relation, 

Y, D (^n)=Q (55) 
n.m 

In case of a generic (even) x (et>en) x (even) lattice with opposite sides 
identified, we give the following sample of the full algebra 

{Cfta,C(n+la)a} = (~ 1)'" 1 si 9 n ( l ) £ ab& e \nb b rn)£ = i' 1 )^ si 9 n ( l )C(n,n+la)a, (56) 

where parentheses denotes symmetrization. 

It is possible to interpret the non-vanishing dirac brackets ([56]) if we treat 
them as constraints. They are then, as we have seen, first class and we may 
write the C(s,m)a' s using (|28| ) as 

C(«, m )a = (4>ij(xft) ~ ^(^KSc^S^a (57) 

we realize that 4>ij(xfi) must be the same for all n's if the C(n jm )a's are to vanish, 

4>ij( x n) = 4>ij( x m) for all n and rfi (58) 

Hence, all vertices are identical, there is only one independent vertex. That is, 
if we consider the " non-closedness of the original algebra" as constraints. 
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4 A four dimensional formulation 



In the past two sections we worked out some details of a Hamiltonian formulation 
of a discretization of general relativity based on distributional fields of the form 
of (2). In this section we would like to sketch out a fully four dimensional 
lagrangian formulation which differs from the one we have just developed in 
that time, as well as space, will be treated discretely. 

To do this we introduce into a four dimensional manifold M., representing 
spacetime, an arbitrary hypercubic lattice. By extending the definitions (2) and 
(15) we will require that the self-dual curvature has support on the faces of 
the lattice, and the self-dual connection has support on the three dimensional 
surfaces of the lattice. We will see that the result of these assumptions is a 
natural reduction of the Lagrangian theory derived from the Capovilla-Dcll- 
Jacobson form (6) to a Lagrangian form of a theory with a countable set of 
degrees of freedom. 

We will follow what we did in the three dimensional case and label the sites 
of the four dimensional lattice by x^, the links by 73a, the faces by <S- a g and 
the three dimensional surfaces by ^ a g e - All of the indices are now four dimen- 
sional and refer to the object that is gotten by moving in the indicated positive 
directions from the site n. Also, the labels of the two and three dimensional 
surfaces are assumed to be antisymmetric in the indices. Assuming a hypercubic 
structure, it is easy to label the faces making up the boundaries of the three 
dimensional surfaces, they are 

nabc nab n+c.ab nac n+o,ac nbc n+a,bc V / 

With these labcllings, we then introduce a Lie algebra valued one form that 
has support on the three dimensional surfaces of the lattice. It is parametrized 
by a Lie algebra element a 8 _ - attached to each three dimensional surface, 1Z~. ;~, 

•> o nabc 1 nabc' 

and it is written as, 

Al(x) = i £ J d 3 ll b n % £ (p)e abcd 6* (x, K n&ie (p)) 4 aSe . (60) 

nabc 

Here, a, b, c, d are four dimensional spacetime indices and p are three coordinates 
on the surface. 

Following the derivation of (18) it is straightforward to show that the curva- 
ture associated with this connection is well defined and distributional, and has 
support on the two dimensional faces of the lattice. It is written as, 

Kb(-) = ± E/ d*S£ b (a)e abcd 5* (x,S n&i (a)) b^, (61) 

nab 

where the Lie algebra elements b\~ associated to each face of the lattice are 

" nab 
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given by 

//. , y (a. , a. . ) ■ -<-> y< , 

nab / 4 \ nabc n—c,abc J / j a>o 



Cy£a,b 



nabc n—c,aoc 



naba n—d,ab 



cd 



(62) 

It is now straightforward to plug (|6l]) into the CDJ form of the action (6) 
and find that the reduced action is 



ncd T n ij 



(63) 



'ibed 



Here q)2 — ^(xfiY 3 is the matrix of scalar fields that satisfies (7) and (8) and 
b\ - is the average of the b\ - 's on the four surfaces in the ab plane that touch 

nab nab 

the site n, 

bL, = ]\bL- h + bi . .,+61 s . s +6l - S J (64) 

nao ^. [ nab n—a,ab n — b,ab n—a — b : ab\ - ' 



5 Conclusions 

In this paper we have introduced two new approximation schemes for general 
relativity which are based on reductions to spaces of solutions based on finite 
dimensional spaces of solutions. Wc believe that it is likely that either the La- 
grangian or the Hamiltonian formulation proposed here could be used to provide 
arbitrary good approximations to solutions to Einstein's equations. However, 
there are a number of things that need to be checked in order to insure that 
either leads to a useful approximation scheme for either classical or quantum 
gravity. We would like to close by listing what we think remains to be done to 
establish the usefulness of the formulations introduced here. 

1) The study of the algebra of the full constraints must be completed. It is 
not necessary that the algebra be first class for a reduction such as the one given 
here to provide a useful approximation method; in the case that the algebra is 
second class the reduction must be understood as constituting a partial gauge 
fixing. However, to go ahead, it is important to know what the algebra is. 

2) The reality conditions must be formulated for the reduced theory, in both 
the Hamiltonian and Lagrangian case. 

3) It may be more convenient, for some purposes, to use a simplicial rather 
than a cubic or hypercubic lattice. In this case it may be useful to introduce a 
dual lattice, whose faces are to be in one to one correspondence with the lines 
of the lattice (in the three dimensional case.) The faces of this dual lattice may 
then be assigned areas which, according to the results described in section 4 
are given by |e/|, where I labels the faces (and the lines of the original lattice.) 
The simplices of this dual lattice can then be taken to constitute a piecewise 
flat Regge manifold whose edge lengths are determined from the areas of the 
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faces. In this way one can reformulate the system described here as a canon- 
ical formulation of the Regge calculus in which the three metric variables are 
treated conventionally, but the conjugate variables are the distributional self- 
dual curvature and connections, which live on a lattice dual to the simplicial 
lattice. 

Furthermore, the Lagrangian formulation for simplicial lattices then appears 
to be a generalization of the four dimensional topological quantum field theories 
of Ooguri]3^| and Crane and Yetter fi4[ . This is a connection that deserves 
further exploration. 

4) It is interesting to speculate the extent to which this formulation could 
be useful for some approach to quantum gravity. It does seems unlikely that 
a direct quantization of the hamiltonian system described in section 2 could 
be useful, given that direct imposition of the diffcomorphism constraints in the 
continuum already reduces quantum gravity to something very much like a fi- 
nite system. On the other hand, it may be that the Lagrangian formulation 
described in section 4 could be the basis for a path integral quantization. Espe- 
cially in the light of the recent impressive progress in four dimensional Monte 
Carlo simulations involving simplicial approximations, it would be interesting to 
investigate the duality suggested by the previous remarks between the present 
formulation and a Regge calculus formulation. It may also be that a direct at- 
tack on the path integral, given the action ( |63| ) may yield interesting results. As 
in any path integral formulation, the key problem that must be solved before 
reliable results can be extracted from the path integral is to find the correct 
measure of integration. 

It is clear that a great deal of work remains to be done to establish whether 
or not the formulations proposed here may be useful for practical calculations in 
classical or quantum gravity. We hope that the results established here suffice 
to justify this further work. 

Appendix. Observables and the meaning of dis- 
tributional geometries 

In this appendix we would like to return to the question of how distributional 
geometries of the kind we have been discussing can be used to approximate arbi- 
trarily well smooth geometries. The idea of using such distributional geometries 
to approximate smooth geometries arose during work on the loop representation 
approach to quantum gravity [^2), when it was realized that the classical limits 
of loop states could be associated with distributional geometries of the form (2). 
It was then discovered that certain classical observables, which are functions of 
the classical fields, could be extended unambigously so that they were defined 
on an extension of the phase space of general relativity that includes distribu- 
tional fields of the form of (2) and (16). Now, this extension cannot be made for 
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all observables, indeed, local observables, such as the metric at a point, cannot 
be given an unambiguous meaning for the distributional geometries. However, 
there are other observables, which are non-local which can be defined in such a 
way that they are meaningful when evaluated on the distributional geometries. 

In references [^2], ^4|, one can find discussions of the quantum operators 
corresponding to these quantities. Very similar considerations apply to the ex- 
tension of the observables to classical distributional geometries, as these results 
have never appeared we describe here the details for one observable, which is 
the area observable. 

The main idea is to use a distributional geometry to approximate a smooth 
geometry by mimicing the weave construction of the quantum theory. For ex- 
ample, if one wants to approximate a metric, q a b 1 which is slowly varying on a 
scale L, up to an accuracy of a << L, one can do the following. We require a set 
of curves a, such that the distributional geometry E%, defined by (1), satisfies 
the following requirement: 

For every surface 5, whose area, given the metric q a h, which we will call 
A[S, q] is larger than L 2 , and whose extrinsic curvature is bounded by l/L 2 , we 
require that 

\A[S,E a }-A[S,q)\<^-^ (65) 

If this is the case then we say that the distributional geometry E^ area-approximates 
the metric q a b on the scale a. 

We may note that the problem of approximating a smooth geometry by a 
distributional one is closely connected to the problem of taking the classical 
limit of a loop state of the quantum theory and showing that it approximates a 
classical metric^, These two problems are closely related because the loop 
states are eigenstates of certain operators that measure the three metric, such 
that the eigenvalues involve the distributional frame field associated to a loop by 
(1). Thus, the problem of extending an observable from smooth to distributional 
geometries is closely related to the problem of constructing good operators for 
those same quantities in the quantum theory through a regularization procedure. 

We would thus like, in the remainder of this section, to show how the area 
operator can be extended to the distributional geometries. From the results it 
will be clear that it is easy to solve the problem of constructing distributional 
geometries that area-approximate any smooth geometry arbitrary well. Analo- 
gous results hold for other observables, including those involving the self-dual 
connections and frame fields, but we do not give them explicitly here. 

Let us consider, then, an arbitrary two dimensional surface in M, which we 
will denote S. The coordinates of S are given by S a (a a ), where the two a a , 
a = 1,2 are coordinates on the surface. Let us the consider the observable A[S], 
which assigns to each surface S its area induced from the three metric on M. 
One would not normally think that this observable could be extended to distri- 
butional geometries. However, we will now show that because our distributional 



20 



frame fields are also densities, A[S] can be defined by a certain procedure, which 
allows it to extend unambiguously to distributional loop geometries. The result 
of this will be that, despite the fact that distributional geometries of the form 
of (2), are defined only on a set of measure zero, the observables A[S] are, when 
nonzero, finite. They assign finite areas to surfaces that cross the loop a. 

We begin by writing the usual expression for A[S] in the case of a smooth, 
nondegenerate geometry, 

A[S] = I Vh, (66) 
Js 

where h is the determinant of the induced two metric, h a b = q a b — n a rib, where 
n a is the unit normal. A simple calculation show that h = q n a n\,. Now, it 
is not hard to show that q cannot itself be extended to distributional loop 
geometries. This can be demonstrated by a direct extension of the argument 

for the non-existence of a unambigous renormalized operator for q in the loop 
representation^^ |2j| As a result, we must construct the area (j6^) through 
a limit that does not need this extension. To do this, let us divide the surface 
up into N disjoint regions Si, such that S — Uj«Si. We then have 

A[S}=J2-A[Si} (67) 

i 

We will proceed by introducing an approximation for the square of «4[<Sj] which 
becomes exact in the limit of infinitesimal surfaces. This is, 

AlpproxlSi] = J d 2 S? b J d 2 S' i cd T**{S,S')ab cd (68) 

where T**(x,y) a f > 0( j = e a be^cdfT e ^ (x , y) . To show that this approximates the 
area of the surface element for small surfaces, we use the facts that in the 
limit T ab (S,S') w <f . We may invert the relation h = q n a n,b to find that 

q = hn a n b — r ab where r ab rib = 0. An infinitesimal element of area is given by 
dA = d 2 S ab Vhn a e abc , from which it follows that, 

dA 2 = d 2 S ab d 2 S cd e ab ee cd fq f (69) 

For smooth fields, this is then equal to (|66| ) in the limit of small areas. We 
may then consider the limit in which we divide the surface up into smaller and 
smaller elements, so that N — > oo. It then follows that, 

N 

A W = lim, E \l A lv P roASil (70) 

For smooth, nondegenerate metrics, this is a long way round to go to define the 
area. But, as we will now show, this particular definition extends to classical 
distributional loop geometries. 
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We then evaluate ( |69| ) for a distributional loop geometry given by (1). We 
find, 

Al pprox [Si] = a 4 f d 2 Sf <f da c (s)5 3 (S,a(s))e abc f d 2 S[ de <f da f (t)S 3 (S' ,a(t))e def 



xTr 

The expression 



w(s)U^ s ,(0,ir)w(t)U^ <s ,(n,2ir) (71) 



J d 2 Sf j> da c 8 3 {S l a)e abc = I[S l7 a] (72) 

is the intersection number of the curve with the surface. It is zero unless they 
intersect, in which case it is equal to ±1 depending on the orientations. Now, 
as we take the limit of N — > oo we will pass a point at which the absolute value 
of the intersection number of a with each surface element is at most one. At 
that point we have 

Al„&] = (a 2 \w( S *)\I[S,a]) 2 (73) 

where s* is the intersection point of the curve with the z'th surface element. At 
this point further subdivisions do not change the value of the sum in (f70|), so 
the limit is equal to 

A[S]=a 2 ^2\^)\ (74) 

i 

where the sum is over the intersection points of the curve with the whole surface. 

Given this result, it is clear how to construct a curve a such that the dis- 
tributional geometry E%, area-approximates a given smooth metric q a b- Indeed, 
the construction mimicks the one that has been already given in the quantum 
case H| M ; the main idea is to arrange the loops so that the sum of the 
contributions in (|70| ) is equal to the area of each surface. 
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